home *** CD-ROM | disk | FTP | other *** search
/ IRIX 6.2 Development Libraries / SGI IRIX 6.2 Development Libraries.iso / dist / complib.idb / usr / share / catman / p_man / cat3 / complib / SPPFA.z / SPPFA
Text File  |  1996-03-14  |  2KB  |  67 lines
  1.  
  2.  
  3.  
  4. SSSSPPPPPPPPFFFFAAAA((((3333FFFF))))                                                            SSSSPPPPPPPPFFFFAAAA((((3333FFFF))))
  5.  
  6.  
  7.  
  8. NNNNAAAAMMMMEEEE
  9.      SPPFA   - SPPFA factors a real symmetric positive definite matrix stored
  10.      in packed form.
  11.  
  12.      SPPFA is usually called by SPPCO, but it can be called directly with a
  13.      saving in time if  RCOND  is not needed.  (Time for SPPCO) = (1 +
  14.      18/N)*(Time for SPPFA) .
  15.  
  16.  
  17. SSSSYYYYNNNNOOOOPPPPSSSSYYYYSSSS
  18.       SUBROUTINE SPPFA(AP,N,INFO)
  19.  
  20. DDDDEEEESSSSCCCCRRRRIIIIPPPPTTTTIIIIOOOONNNN
  21.      On Entry
  22.  
  23.      AAAAPPPP REAL (N*(N+1)/2)
  24.         the packed form of a symmetric matrix  A .  The
  25.         columns of the upper triangle are stored sequentially
  26.         in a one-dimensional array of length  N*(N+1)/2 .
  27.         See comments below for details.
  28.  
  29.      NNNN INTEGER
  30.         the order of the matrix  A .  On Return
  31.  
  32.      AAAAPPPP an upper triangular matrix  R , stored in packed
  33.         form, so that  A = TRANS(R)*R .
  34.  
  35.      IIIINNNNFFFFOOOO INTEGER
  36.         = 0  for normal return.
  37.         = K  if the leading minor of order  K  is not
  38.         positive definite.  Packed Storage The following program segment will
  39.      pack the upper triangle of a symmetric matrix.
  40.         K = 0
  41.         DO 20 J = 1, N
  42.         DO 10 I = 1, J
  43.         K = K + 1
  44.         AP(K) = A(I,J)
  45.         10    CONTINUE
  46.         20 CONTINUE LINPACK.  This version dated 08/14/78 .  Cleve Moler,
  47.      University of New Mexico, Argonne National Lab.  Subroutines and
  48.      Functions BLAS SDOT Fortran SQRT
  49.  
  50.  
  51.  
  52.  
  53.  
  54.  
  55.  
  56.  
  57.  
  58.  
  59.  
  60.  
  61.  
  62.  
  63.                                                                         PPPPaaaaggggeeee 1111
  64.  
  65.  
  66.  
  67.